7-Z box corrections to weak charges of heavy nuclei in atomic parity violation 
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We present a new dispersive formulation of the 7Z box radiative corrections to weak charges of 
bound protons and neutrons in atomic parity violation (APV) measurements on heavy nuclei such as 
133 Cs and 213 Ra. We evaluate for the first time a small but important additional correction arising 
from Pauli blocking of nucleons in a heavy nucleus. Overall, we find a significant shift in the 7Z 
correction to the weak charge of 133 Cs, approximately 4 times larger than the current uncertainty 
on the value of sin 2 Ow, but with a reduced error compared to earlier estimates. 
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In the search for physics beyond the Standard Model 
one of the most important indirect methods involves a 
high precision test of the evolution of the Weinberg an- 
gle with scale. In particular, the comparison between the 
value of sin 2 6w measured at the Z pole at LEP and the 
value extracted from parity violation in atomic systems 
has provided a very strong confirmation of the radiative 
corrections calculated within the Standard Model. The 
atomic system with the most accurate current measure- 
ment is the parity-violating S—S transition in neutral 
The predicted enhancement factor of the 
APV effect in the S-D transition in 213 Ra + is about 50 
times larger Q . It is critical that the calculations of the 
radiative corrections for these systems incorporate the 
latest theoretical developments and match the precision 
of the experimental data. 

Driven by the demand for accurate radiative correc- 
tions for high energy parity- violating electron scattering, 
notably the Q W eak experiment at Jefferson Lab 0] , there 
have recently been new evaluations of the 7^ box dia- 
gram 5 9]. In particular, Blunden et al. Q developed 
a formulation of this correction in terms of dispersion 
relations and moments of the inclusive 7Z interference 
structure functions. This approach provides a systematic 
method for improving the accuracy of the calculation. 

In this Letter we use the dispersive relations methods 
to compute a new value for the radiative correction as- 
sociated with the 7Z box diagram for both protons and 
neutrons in 133 Cs and 213 Ra. Compared with previous 
estimates that were computed some 3 decades ago fioj ]. 
the new corrections give jZ contributions that are 15% 
smaller for free nucleons, and 20% smaller for nucleons 
bound in a heavy nucleus (due to Pauli blocking). This 
shifts the theoretical value of Qw(Cs) from —73.14(6) to 
—73.26(4), which agrees well with the current experimen- 
tal value of -73.20(35). 

Including radiative corrections, the weak charges of the 
proton and neutron can be written as 
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with s 2 = sin 2 9 W (M 2 Z ) = 0.23116(13) in the MS scheme. 
A e and A'„ are Zee and jee vertex corrections, given in 
Refs. [hHHI, and expressions for the universal parame- 
ters p and k(Q 2 ) are given in Ref. [12 1. 

The proton weak charge Qiy is sensitive to the running 
of the Weinberg angle, as embodied in k(Q 2 ), whereas the 
neutron weak charge QIV 1S primarily sensitive to p. At 
the one-loop level, p has a quadratic dependence on the 
top quark mass. This dependence is modified by signifi- 
cant higher-order QCD and electroweak corrections [IH , 
some of which have been evaluated to four loops. 

The correction k(Q 2 ) includes boson self-energy contri- 
butions from "fZ mixing. In particular, it has a hadronic 
uncertainty from the quark contributions to fermion 
loops, denoted by Ak^ for five quark flavors, which 
is strongly correlated with the analogous contribution 



A^had to tne running of a(Q 2 ). A reduction of Aa 



(5) 
had 



to 



0.02772 in the most recent analysis implies a corrcspond- 



ing reduction in Aft^ ad compared with Ref. [Ill ], leading 
to Ak^s 2 = 7.87(8) x 10" 3 . Using the most recent Stan- 
dard Model parameters [Hj|, we obtain p = 1.0006(2) 
and k(0)s 2 = 0.23807(15), in essential agreement with 
the values quoted in Ref. [l4j. 

The remaining terms in Eqs. ([T]) and ([2]) arise from the 
(numerically dominant) WW, ZZ, and 7Z boxes. We 
take the expressions for WW and ZZ boxes from Erlcr 
et al. [ll| , including the leading order QCD correction to 
the original expressions of Ref. [Tcj . The box diagrams 
with two heavy bosons are dominated by high momentum 
scales By contrast, the 7Z box diagram contains 

both high and low momentum scales, and is therefore 
sensitive to hadronic corrections. Following convention 
we write the jZ contribution for protons and neutrons 
in terms of a parameter B N , 
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= 5/3(4/3) for the proton 
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where v e 

(neutron). A free-quark model gives B p = \n(M^/m 2 ) 
3/2, where m is an undetermined hadronic mass scale 
(such as a constituent quark mass); however, this does 
not adequately describe either the long-range or short- 
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TABLE I. Weak charges of the proton, neutron, 133 Cs, and 
213 Ra, comparing the previous MS estimates OH EH with our 
new results for free and bound nucleons. 





MS 


free nucleons 


bound nucleons 


B p 


11.8(1.0) 


9.95(40) 


9.36(40) 




11.5(1.0) 


9.82(40) 


9.32(40) 


Qw 


0.0716(8) 


0.0708(6) 


0.0705(6) 


Qw 


-0.9882(4) 


-0.9888(2) 


-0.9890(2) 


Qw (Cs) 


-73.14(6) 


-73.23(4) 


-73.26(4) 


Qw(Ra.) 


-117.22(10) 


-117.37(7) 


-117.42(7) 



with (M + m^) 2 < W 2 < 4 GcV 2 ; and (iii) deep-inelastic 
scattering (DIS), with W 2 > 4 GeV 2 . Contributions 
from region (i) depend on the nucleon's elastic magnetic 



G^j and axial- vector G A 



Z,N 



form factors, 



F N{el) (x 



-GZ(Q 2 )G% N (Q 2 )x5(l-x). (7) 



We set GZ(Q 2 ) = p N F v {Q 2 ), and G Z A ' N (Q 2 ) = 
—g A FA(Q 2 ), with fi N the nucleon magnetic moment, 
and g p A = —g A = 1.267. A dipole Q 2 dependence, 
F v ,a(Q 2 ) = + Q 2 / a v,a) 2 , suffices for both V and A 
form factors, with Ay = 6.84 GcV, and A^ = 1.0 GcV. 
More sophisticated form factors give essentially identical 
numerical results. 

For the resonance contributions we use the 
parametrizations of the transition form factors from 
Lalakulich et al. [l6j |. but with modified isospin factors 
appropriate to jZ. These form factors have been fit to 
pion production data in v and v scattering, and include 
the lowest four spin-1/2 and 3/2 states. 

For the DIS region, we divide the Q 2 integral of Eq. (J6]) 
The error in Q p w (excluding yZ boxes) arises from k(0)s 2 i nto a low-Q 2 part Q 2 < Q§, where the structure function 
(±0.0006) and the WW boxes (±0.0001), while the cr- t?n ; „ Q i Q +;,, Q i, r ,„-,w.„„ rrl „„a „ v,^v. n2 tr& 
tot in QVrr arises from p (±0.0002). The results in the 
column labeled MS in Table U use the most recent esti- 
mates of B p and B n from Refs. [IHEij]. For Cesium, we 
have Qw (Cs) — 55Q^±78QiV- Adding the independent 
errors in quadrature (the errors in B p and B n are not in- 
dependent), we find —73.14(6). The yZ boxes contribute 
±0.052 to this total, while all other errors combined con- 
tribute ±0.037. 



range behavior of the yZ box. Marciano and Sirlin (MS) 
[lfjj give a more refined estimate by separately modeling 
the low and high energy contributions, but an equivalent 
logarithmic dependence on the hadronic mass parameter 
m remains. 

Using Eqs. ([T]) and ([2]), we find numerically 

Q p w = 0.0664(6) + 0.00044 B p , (4) 
Q" w = -0.9922(2) + 0.00035 B n . (5) 



Fj^ z is relatively unknown, and a high-Q 2 part (Q 2 > 
Qg), where at leading order the structure function can 
be expressed in terms of valence quark distributions [l4T ] , 



f 3 n T ] ^Q 2 ) 
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To proceed beyond the work of MS 10( we use the 
dispersion methods developed in Refs. [E-Q, writing the 
imaginary part of n lZ [E) in terms of structure func- 
tions Ff z (i = 1, 2, 3) that can be obtained from inclu- 
sive lepton-nucleon scattering. A dispersion integral over 
energy then gives the real part of a jZ (E), which con- 
tributes to the weak charge. The energy dependence of 

"~[. At E = 0, relevant 
for APV, only the axial-vector ZN coupling involving 

3,-fZ 



At high Q 2 the yZ contribution can be expanded in pow- 
ers of x 2 /Q 2 , yielding a series whose coefficients are struc- 
ture function moments of increasing rank, 

±r d0 ._*(Q 2 )v e (Q 2 ^ 
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□ 7 z(-E) was evaluated in Refs. 
PV 

the F^ lZ structure function gives a nonzero contribu- 
tion. Performing the dispersion energy integral analyti- 
cally, the real part of O jZ can be written 

2 r ,„o a(Q 2 )v e (Q 2 ) 



where the n-th moment of the F7 structure function is 
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where 7 = y/l + 4M 2 x 2 /Q 2 and x = Q 2 /(W 2 - M 2 + 
Q 2 ), with W the invariant mass of the intermediate 
hadronic state. Following Ref. [§] , we include the running 
with Q 2 of a{Q 2 ) and v e (Q 2 ) = 1 - 4k(Q 2 ) s 2 in Eq. © 
due to boson self-energy contributions. Both quantities 
vary significantly over the relevant Q 2 range. 

The contributions to O jZ can be split into three kine- 
matic regions: (i) elastic, with W 2 — M 2 ; (ii) resonances, 



Numerically, the n = 1 moment dominates, with the n > 
3 contributions to the integral of Eq. © less than 0.1%. 

The lowest moment is in fact the 7Z analog of the 
GLS sum rule for vN DIS, which at leading order 
counts the number of valence quarks in the nucleon. The 
corresponding quantity for 7 Z is Q p — ^2 q 2e q g q A = 5/3 
for the proton and Q n = 4/3 for the neutron. Including 
the next-to- leading order strong interaction correction in 
the MS scheme, the n = 1 contribution is 



M^\Q 2 ) = Q N 



1 - 



a s (Q 2 ) 
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Combined with Eq . this is identical to the high en- 
ergy result of MS [10j, but with Q replacing the arbi- 
trary hadronic mass parameter m. In our case the scale 



3 



Qo corresponds to the momentum above which a partonic 
representation of the nonresonant structure functions is 
valid, and above which the Q 2 evolution of parton distri- 
bution functions (PDFs) via the Q 2 evolution equations 
is applicable. We vary Qo between 1 and 2 GeV 2 , which 
coincides with the typical lower limit of recent sets of 
PDFs. 

The contribution for Q 2 < Q 2 can in principle be ob- 
tained from data. There is limited information on F^.w 
from neutrino scattering, but little or no existing data on 
F3 n z at low Q 2 . As in Ref. @, we use two different mod- 
els to smoothly interpolate in Q 2 between Qq and 0: one 
vanishes in the Q 2 — > limit, and the other approaches a 
constant. The differences between the models are of the 
order 5-15%. 

The relation between proton and neutron contributions 
from the different kinematic regions is 
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The near equality of the resonance contributions (within 
3%) is due to the dominance of isovector resonances, 
which contribute equally for protons and neutrons. 

The full results are summarized in Table HU The largest 
contributions come from the DIS region. Fortunately, 
the results show only a mild sensitivity to the parame- 
ter Qq that separates the model dependent low-Q 2 ex- 
trapolation from the high-Q 2 partonic region. The main 
part of the uncertainty arises from the model dependence 
of this extrapolation to low Q 2 . We therefore assign 
B p = 9.95(40), equal to the average of the four values 
in Table [IIJ with a very conservative error given by the 
DIS contribution from the region Q 2 < 1 GeV 2 . This is 
the value appearing in the column labeled 'free nucleons' 
of Table [IJ together with the neutron contribution using 
Eq. (Ull). 

Because the nucleons in a heavy nucleus are bound, 
properties such as their weak charge can differ from those 
of free nucleons. In particular, for the elastic "fZ box 
contribution, transitions to occupied states are forbidden 
by the Pauli exclusion principle. To estimate this effect 



TABLE II. Contributions to B p from different kinematic re- 
gions for two different values of the matching scale Qq. The 
range of values for the DIS (Q 2 < Qq) contribution is for the 
two models described in Ref. 





1 GeV 2 


2 GeV 2 


elastic 


1.47 


1.47 


resonance 


0.59 


0.59 


DIS 






Q 2 > Qo 


7.50 


7.05 


Q 2 < Qo 


0.42-0.48 


0.78-0.82 


Total 


9.98-10.04 


9.89-9.93 



of the nuclear medium on the weak charges, we consider 
the expression 



= -Vet* N 9% r 'dQ 2 F v (Q 2 )F A (Q 2 )f(Q 2 ) 
n Jo 
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where 
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Here v e is taken at an appropriate low-momentum scale, 
and the Q 2 dependence of the Z propagator has been 
dropped. 

Pauli blocking is important because the integrand in 
Eq. (fTB")) is heavily weighted towards low Q 2 , 
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Since the form factors introduce corrections of order 
Q 2 /A 2 , the dominant low-Q 2 contribution is largely in- 
dependent of nuclcon structure. 

To allow for Pauli blocking of the intermediate nuclcon 
in a heavy nucleus we use the Fermi gas model, where 
nucleon states are occupied below the Fermi momentum 
kp (typically ~ 0.25 GeV). This estimate should suf- 
fice in a heavy nucleus like 133 Cs with a low surface to 
volume ratio. For momentum transfer q to a bound nu- 
cleon of momentum p, we must exclude from the integral 
of Eq. (|13j) all intermediate nucleon states of momen- 
tum \p + q\ < kp. Introducing the occupation number 
n p = 0(kp — p), we therefore have a factor 



C(q) 



J d 3 pi 
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to be folded into the integrand of Eq. (fTB")) , with 
J d 3 pn p = (47r/3)/cp. This represents the fractional vol- 
ume of occupied states that cannot reach the Fermi sur- 
face for a given value of q. Since kp -C M 2 , the affected 
states have nonrelativistic energies, and so Q 2 ps q 2 . 
From simple geometry, we find 



C{q) 



1 2 ( 3 2 k r 
0. 



< q < 2k F , 

q > 2k F . 



(17) 



The expression 1 — C(q) is also the Coulomb sum rule 
for longitudinal quasielastic scattering in a nonrelativistic 
Fermi gas. 

The Pauli blocking effect can be introduced as a cor- 
rection 1 — A(/cp) to the value of in Eq. (|13|) . with 



J 2kF dQ 2 F V (Q 2 )F A (Q 2 ) f(Q 2 )C(Q) 
J °° dQ 2 F V {Q 2 )F A (Q 2 ) f(Q 2 ) 
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representing the fractional contribution from the ex- 
cluded states. This correction depends on kp, and has 



4 



TABLE III. The parameters of the proton and neutron 
Woods-Saxon densities p(r) for 133 Cs, together with the mean 
square radii. The Fermi momenta kp(0) are determined from 
the central nucleon densities, with (fc^) the corresponding 
mean. 



c (fm) 


a (fm) 


(r 2 > (fm 2 ) 


k F (0) (GeV) 


(k F ) (GeV) 


p 5.8895 


0.4010 


23.03 


0.241 


0.218 


n 5.9482 


0.4946 


24.61 


0.266 


0.236 



only a weak dependence on nucleon structure through the 
form factor parameters A, with other parameters drop- 
ping out in the ratio. In the numerator, the leading terms 
are of order O(kp) and 0{kp) 1 with the form factor cor- 
rections only appearing at order 0{k F ) and higher. 

The values of kp for protons and neutrons are taken 
from Hartree-Fock calculations for 133 Cs, which repro- 
duce the experimental charge density [18|. Fitting the 
Hartrcc-Fock proton and neutron distributions to a stan- 
dard Woods-Saxon form p(r) = po/(l + cxp[(r — c)/a]), 
normalized to Z and N, respectively, leads to the pa- 
rameters given in Table IIIII From these one can compute 
the Fermi momentum in the local density approximation, 
p = kp/(3% ). The central and average values of kp are 
given in Table Hill The latter, which are used in our cal- 
culations, are consistent with those obtained in fits to 
experimental quasielastic electron scattering data using 
a simple Fermi gas model 

We find the Pauli blocking correction factor 1 — A(/cp) 
is approximately linear in kp over the range 0.2—0.3 GeV, 
and well- approximated by the expression 



l-A(kp) « 0.83- 1.04 kp, 



(19) 



with kp in GeV. Specifically, using the volume-averaged 
values of (kp) in Table Hill we find a correction factor of 
0.61 to and 0.59 to • For reasonable values 

of kp up to 0.27 GeV in a very heavy nucleus, the Pauli 
blocking factors will therefore fall into the narrow range 
0.55 — 0.60, suggesting a relatively insignificant variation 
in the total value of n^ z for nuclei beyond 133 Cs. 

The effect on Qw(Cs) is shown in final column of Ta- 
bled] There is a small shift of —0.03 compared to the free 
nucleon values, giving —73.26(4). The individual terms 
contributing to the uncertainty are: k(0)s 2 (±0.033), 
WW boxes (±0.006), p (±0.016), and 7 Z (±0.021). Also 
included in Table U is our theoretical value for the weak 
charge of 213 Ra. 

In summary, we have computed the effect of ex- 
change corrections on the weak charge of heavy nuclei 
such as 133 Cs and 213 Ra, using a recently developed for- 
malism based on dispersion relations and an expansion 
of the ^Z interference structure function moments. The 



results improve earlier estimates based on a quark model 
description of the jZ box contributions, allowing a signif- 
icant reduction in the theoretical uncertainty. Compared 
with the pioneering early estimates of Marciano and Sir- 
lin [roj . the new corrections enhance the magnitude of 
the weak charge by w 0.16%, which is approximately 4 
times larger than the current uncertainty on sin 2 8\y , and 
will therefore affect future high-precision determinations 
of the weak angle. 
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